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Abstract 

Using the coupled cluster method (CCM) we study the zero-temperature phase diagram of a 
spin-half Heisenberg antiferromagnet (HAF), the so-called J1-J2 model, defined on an anisotropic 
two-dimensional lattice. With respect to an underlying square-lattice geometry the model contains 
antiferromagnetic (Ji > 0) bonds between nearest neighbors and competing (Jg > 0) bonds be- 
tween next-nearest neighbors across only one of the diagonals of each square plaquette, the same 
diagonal in every square. Considered on an equivalent triangular-lattice geometry the model may 
be regarded as having two sorts of nearest-neighbor bonds, with J2 = kJi bonds along parallel 
chains and Ji bonds providing an interchain coupling. Each triangular plaquette thus contains 
two Ji bonds and one J2 bond. Hence, the model interpolates between a spin-half HAF on the 
square lattice at one extreme {k = 0) and a set of decoupled spin-half chains at the other {k — > 00), 
with the spin-half HAF on the triangular lattice in between at k = 1. We use a Neel state, a 
helical state, and a collinear stripe-ordered state as separate starting model states for the CCM 
calculations that we carry out to high orders of approximation (up to LSUB8). The interplay 
between quantum fluctuations, magnetic frustration, and varying dimensionality leads to an inter- 
esting quantum phase diagram. We find strong evidence that quantum fiuctuations favor a weakly 
first-order or possibly second-order transition from Neel order to a helical state at a first critical 
point at = 0.80 lb 0.01, by contrast with the corresponding second-order transition between the 
equivalent classical states at Kd = 0.5. We also find strong evidence for a second critical point at 
= 1-8 ± 0.4 where a first-order transition occurs, this time from the helical phase to a collinear 
stripe-ordered phase. This latter result provides quantitative verification of a recent qualitative 
prediction of Starykh and Balents [Phys. Rev. Lett. 98, 077205 (2007)] based on a renormalization 
group analysis of the J1-J2 model that did not, however, evaluate the corresponding critical point. 

PACS numbers: 75.10.Jm, 75.30.Gw, 75.40.-s, 75.50.Ee 
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I. INTRODUCTION 

Two-dimensional (2D), spin-1/2, Heisenberg antiferromagnets (HAFs) have been much 
studied in recent years. The interplay between (either dynamic or geometric) frustration and 
quantum fluctuations in determining the ground-state (gs) phase diagram of such models 
has been of particular interest. While such models are well understood in the absence 
of frustration/ this is not the case for frustrated systems, for which the zero-temperature 
(T = 0) phase transitions between magnetically ordered quasiclassical phases and novel 
(magnetically disordered) quantum paramagnetic phases^i^ have become the subject of great 
recent interest. A particularly well studied such model is the frustrated J1-J2 model on 
the square lattice with nearest-neighbor (NN) bonds (Ji) and next-nearest-neighbor (NNN) 
bonds {J2), for which it is now well accepted that there exist two phases exhibiting magnetic 
long-range order (LRO) at small and at large values of a = J2/J1 respectively, separated 
by an intermediate quantum paramagnetic phase without magnetic LRO in the parameter 
regime < a < ac2, where ~ 0.4 and ~ 0.6. For a < the gs phase exhibits Neel 
magnetic LRO, whereas for a > a^^ it exhibits coUinear stripe LRO. We have recently studied 
this 2D spin-1/2 model exhaustively by extending it to include anisotropic interactions in 
either real (crystal lattice) spaced or in spin space.- We showed in particular how the coupled 
cluster method (COM) provided for this highly frustrated model what is perhaps now the 
most accurate microscopic description. The interested reader is referred to Refs. and 
references cited therein for further details of the model and the method. 



II. THE MODEL 

In the light of the above successes we now apply the CCM to the seemingly similar 2D 
spin-1/2 J1-J2 niodel that has been studied recently by other means .-'^>^'^>^>^ Its Hamilto- 
nian is written as 

i7 = Ji ^ Si ■ + J2 ^ Si • Sfc (1) 

where the operators Si = (sf , s^, sf) are the spin operators on lattice site i with sf = s(s + 1) 
and s = 1/2. On the square lattice the sum over (i, j) runs over all distinct NN bonds, 
but the sum over [i, fc] runs only over one half of the distinct NNN bonds with equivalent 
bonds chosen in each square plaquette, as shown explicitly in Fig. [H (By contrast, the J\- 
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FIG. 1: (Color online) Ji-J'2 model; — Ji; Jg; (a) Neel state, (b) spiral state, (c) stripe state. 

J2 model discussed above includes all of the diagonal NNN bonds.) We shall be interested 
here only in the case of competing (or frustrating) antiferromagnetic bonds Ji > and 
J2 > 0, and henceforth for all of the results shown we set Ji = 1. Clearly, the model 
may be described equivalently as a Heisenberg model on an anisotropic triangular lattice 
in which each triangular plaquette contains two NN Ji bonds and one NN bond. The 
model thus interpolates continuously between HAFs on a square lattice {J 2 = 0) and on a 
triangular lattice {J 2 = Ji)- Similarly, when Ji = (or J2 ^ 00 in our normalization with 
Ji = 1) the model reduces to uncoupled ID chains (along the chosen diagonals on the square 
lattice). The case ^ 1 thus corresponds to weakly coupled ID chains, and hence also 
interpolates between ID and 2D. We note in this context that the CCM has also been very 
successfully applied to other spin-1/2 HAF models that continuously interpolate between 
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(a) the triangular and kagome lattices;J^ (b) the square and honeycomb lattices; 
and 2D cases;- and (d) 2D and 3D cases.^^ As well as the obvious theoretical richness of the 
model, there is also experimental interest since it also well describes such quasi-2D materials 
as BEDT-TTF crystalsi^ with J^/Ji ^ 0.34-1, and CssCuCl^i^ with J^/Ji ^ 6. 

The J1-J2 model has only two classical gs phases (corresponding to the case where the spin 
quantum number s 00). For J2 < \Ji the gs phase is Neel ordered, as shown in Fig. [T](a), 
whereas for J'2 > |Ji it has spiral order, as shown in Fig. [l](b), wherein the spin direction at 
lattice site (i,j) points at an angle = + {i + j)aci, with ad = cos~^{—jjr) = vr — 0ci- 



The pitch angle (pci = cos' 



-i( Ji 
y2j' 



thus measures the deviaton from Neel order, and it varies 



from zero for 2 Jg/ Ji < 1 to |7r as J2/ Ji 00, as shown in Fig. [3l When = Ji we regain 
the classical 3-sublattice ordering on the triangular lattice with = Ivr. The classical 
phase transition at J2 = ^Ji is of continuous (second-order) type, with the gs energy and 
its derivative both continuous. 

In the limit of large jy Ji the above classical limit represents a set of decoupled ID HAF 
chains (along the diagonals of the square lattice) with a relative spin orientation between 
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neighboring chains that approaches 90°. In fact, of course, there is complete degeneracy at 
the classical level in this limit between all states for which the relative ordering directions 
of spins on different HAF chains are arbitrary. Clearly the exact spin- 1/2 limit should also 
be a set of decoupled HAF chains as given by the exact Bethe ansatz solution. However, 
one might expect that this degeneracy could be lifted by quantum fluctuations by the well- 
known phenomenon of order by disorderM Just such a phase is known to exist in the J1-J2 
model^i^ for values of J2/J1 ^ 0.6, where it is the so-called coUinear stripe phase in which, 
on the square lattice, spins along (say) the rows in Fig. [1] order ferromagnetically while 
spins along the columns and diagonals order antiferromagnetically, as shown in Fig. [H^c). 
We note, however, that a corresponding order by disorder phenomenon, if it exists for the 
present J1-J2 ^lodel, would be more subtle than for its textbook J1-J2 model counterpart, 
as we explain more fully in Sec. V. 

In a recent paper Starykh and Balents^ have given a renormalization group (RG) analysis 
of the spin-1/2 J1-J2 model considered here to predict that precisely such a coUinear stripe 
phase also exists in this case for values of jy Ji above some critical value (which they do 
not calculate). One of the aims of the present paper is to give a fully microscopic analysis of 
this model in order to map out its T = phase diagram, including the positions and orders 
of any quantum phase transitions that emerge. 



III. THE COUPLED CLUSTER METHOD 



The CCM (see, e.g., Refs. 
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2l| and references cited therein) that we employ here 



is one of the most powerful and most versatile modern techniques in quantum many-body 
leory. It has been applied very successfully to various quantum magnets (see Refs. 



theory. It ha 



231] and references cited therein). The method is particularly appropriate 



for studying frustrated systems, for which the main alternative methods are often only of 
limited usefulness. For example, quantum Monte Carlo techniques are particularly plagued 
by the sign problem for such systems, and the exact diagonalization method is restricted in 
practice, particularly for s > 1/2, to such small lattices that it is often insensitive to the 
details of any subtle phase order present. 

The method of applying the CCM to quantum magnets has been described many times 



elsewhere (see, e.g., Refs. [12| 



mm 



2l| and references cited therein) . It relies on build- 
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ing multispin correlations on top of a chosen gs model state |$) in a systematic hierarchy 
of LSUBra approximations for the correlation operators 5* and S that exactly parametrize 
the exact gs ket and bra wave functions of the system respectively as |^) = e'^l'l') and 
= ($|S'e^'^. In the present case we use three different choices for the model state |$), 
namely either of the classical Neel and spiral states, as well as the collinear stripe state. 
Note that for the helical phase we perform calculations for arbitrary pitch angle a = vr — 0, 
and then minimize the corresponding LSUBn approximation for the energy with respect 
to 0, -ELSUBn(0) min cj) = 0LSUBn- Generally (for n > 2) the minimization must be 
carried out computationally in an iterative procedure, and for the highest values of n that 
we use here the use of supercomputing resources was essential. Results for 0LSUBn will be 
given later (Fig. [3]). We choose local spin coordinates on each site in each case so that all 
spins in |$), whatever the choice, point in the negative z-direction (i.e., downwards). 

Then, in the LSUBn approximation all possible multi-spin-fiip correlations over different 
locales on the lattice defined by n or fewer contiguous lattice sites are retained. Clearly, 
in the present case we have a choice whether to consider the model to be defined on the 
square lattice (shown in Fig. [T]) or to consider it on the (topologically equivalent) triangular 
lattice, as discussed in Sec. [Tll Although these two viewpoints are completely equivalent for a 
description of the model, they differ for the purposes of defining the LSUBn approximations. 
Thus, for example, each pair of sites joined by a J2 bond are NNN pairs on the square lattice 
but are NN pairs on the triangular lattice. Hence, such a (NN on the triangular lattice) 
double spin-flip configuration is contained in LSUBn approximations on the square lattice 
only for n > 3, whereas it is contained at the LSUBn level on the triangular lattice for 
n > 2. Whereas both LSUBn hierachies agree in the n —* 00 limit they will differ for finite 
values of n. In general there are clearly more multi-spin-flip conflgurations retained at a 
given LSUBn level on the triangular lattice than on the square lattice, and in the present 
paper we consider only the triangular case. 

The numbers of such distinct fundamental conflgurations on the triangular lattice (viz., 
those that are distinct under the space and point group symmetries of both the Hamiltonian 
and the model state |$)) that are retained for the collinear stripe and spiral states of the 
current model in various LSUBn approximations are shown in Table [H The coupled sets 
of equations for these corresponding numbers of coefficients in the operators S and S are 
derived using computer algebrar^ and then solved^ using parallel computing. We note that 
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TABLE I: Number of fundamental LSUBn configurations ({j f.c.) for the stripe and spiral states of 
the spin-1/2 model, using the triangular lattice geometry. 

tt f.c. 

Method 

stripe spiral 

LSUB2 2 3 

LSUB3 4 14 

LSUB4 27 67 

LSUB5 95 370 

LSUB6 519 2133 

LSUB7 2617 12878 

LSUB8 15337 79408 

such CCM calculations using up to about 10^ configurations or so have been previously 
carried out many times using the CCCM code^^ and heavy parallelization. A significant 
extra computational burden arises here for the helical state due to the need to optimize the 
quantum pitch angle at each LSUBn level of approximation as described above. Further- 
more, for many model states the quantum number = Xlili tie used to restrict 
the numbers of fundamental multi-spin-flip configurations to those clusters that preserve 

= 0. However, for the spiral model state that symmetry is absent, which largely explains 
the significantly greater number of fundamental configurations for the spiral state than for 
the stripe state at a given LSUBn order. Hence, the maximum LSUBn, level that we can 
reach here, even with massive parallelization and the use of supercomputing resources, is 
LSUB8. For example, to obtain a single data point (i.e., for a given value of J2, with Ji = 1) 
for the spiral phase at the LSUB8 level typically required about 0.3 h computing time using 
600 processors simultaneously. 

At each level of approximation we may then calculate a corresponding estimate of the 
gs expectation value of any physical observable such as the energy E and the magnetic 
order parameter, M = — (\l/|sf|^), defined in the local, rotated spin axes, and which thus 
represents the on-site magnetization. Note that M is just the usual sublattice magnetization 
for the case of the Neel state as the CCM model state, for example. More generally it is 
just the on-site magnetization. 
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It is important to note that we never need to perform any finite-size scaling, since all 
CCM approximations are automatically performed from the outset in the infinite-lattice 
limit, — oo, where N is the number of lattice sites. However, we do need as a last step 
to extrapolate to the n ^ oo limit in the LSUBn truncation index n. We use here the 
well-testecy^ii^ empirical scaling laws 

E/N = ao + aiTi'^ + a2n~'^ , (2) 

M = bo + hn-^ + b2n-^ , (3) 

that have given good results previously, for example, for the interpolating triangle-kagome 
HAF^^ and the interpolating square-honeycomb HAF.— We comment further on the accu- 
racy of the extrapolations in Sec. |V] where we present a discussion of our results. 

IV. RESULTS 

We report here on CCM calculations for the present spin-1/2 J1-J2 Kiodel Hamiltonian 
of Eq. ([1]) for given parameters (Ji = 1, Jg), based respectively on the Neel, spiral and 
stripe states as CCM model states. Our computational power is such that we can perform 
LSUBn calculations for each model state with n < 8. We note that, as has been well 
documented in the past,— the LSUBn data for both the gs energy per spin E/N and the on- 
site magnetization M converge differently for the even-ra sequence and the odd-n sequence, 
similar to what is frequently observed in perturbation theory. Since, as a general rule, it 
is desirable to have at least (n + 1) data points to fit to any fitting formula that contains n 
unknown parameters, we prefer to have at least 4 results to fit to Eqs. ([2]) and ([2]). Hence, for 
most of our extrapolated results below we use the even LSUBn sequence with n = {2, 4, 6, 8}. 

We report first on results obtained using the spiral model state. While classically we 
have a second-order phase transition from Neel order (for k, < Kci) to helical order (for 
K > Kc\), where n = J2/J1, at a value Hci = 0.5, using the CCM we find strong indications 
of a shift of this critical point to a value k^^ ^ 0.80 in the spin-1/2 quantum case. Thus, 
for example, curves such as those shown in Fig. [2] show that the Neel model state (0 = 0) 
gives the minimum gs energy for all values of k, < where is also dependent on the 
level of LSUBn approximation, as we also see below in Fig. [31 By contrast, for n > 
the minimum in the energy is found to occur at a value 7^ 0. If we consider the pitch 
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FIG. 2: (Color online) Ground-state energy per spin of the spin-1/2 J1-J2 Hamiltonian of Eq. ([1]) 
with Ji = 1, using the LSUB6 approximation of the CCM with the spiral model state, versus the 
spiral angle 0, for some illustrative values of J2 in the range < J2 < 1.4. For < 0.788 the 
minimum is at = (Neel order), whereas for JI2 > 0.788 the minimum occurs aX (j) = 4>hS\]B& 7^ 0, 
indicating a phase transition at J2 ~ 0.788 in this approximation. 

angle itself as an order parameter (i.e., = for Neel order and 7^ for spiral order) a 
typical scenario for a first-order phase transition would be the appearance of a two-minimum 
structure for the ground-state energy as a function of 0, exactly as shown in Fig. [2] for the 
LSUB6 approximation. Very similar curves occur for other LSUBn approximations. If we 
therefore admit such a scenario, in the typical case one would expect various special points 
in the transition region, namely the phase transition point itself where the two minima 
have equal depth, plus one or two instability points /tj^ and where one or other of the 
minima (at = and 7^ respectively) disappears. In the present case, it is interesting to 
note that the two points and either coincide exactly or are indistinguishable within 
the accuracy of our calculations, thereby indicating that the transition at is rather subtle, 
and perhaps even of second-order rather than first-order type. A close inspection of curves 
such as those shown in Fig. [2] for the LSUB6 case shows that what happens at this level 
of approximation is that for n < 0.788 the only minimum in the ground-state energy is 
at = (Neel order). As this value is approached asymptotically from below the LSUB6 
energy curves become extremely flat near 7^ 0, indicating the disappearance at = of 
the second derivative (PE/dcfP' (and possibly also of one or more of the higher derivatives 
(T'E/dcf)"' with n > 3), as well as of the first derivative dE/d(f). Then, for all values n > 0.788 




9 




FIG. 3: (Color online) The angle (pLSVBn that minimizes the energy i?LSUBra(</') of the spin-1/2 
J1-J2 Hamiltonian Eq. ([1]) with Ji = 1, in the LSUBn approximations with n = {2,4,6,8}, using 
the spiral model state, versus J2. The corresponding classical result (pci is shown for comparison. 
We find in the LSUBn quantum case with n > 2 a seemingly first-order phase transition (e.g., for 
LSUB8 at J2 ~ 0.796 where (/>lsub8 jumps abruptly from zero to about 0.147r), although also see 
the text for a broader discussion of the nature of the quantum phase transition. By contrast, in 
the classical case there is a second-order phase transition at J2 = 0.5. 



the LSUB6 curves develop a secondary minimum at a value 7^ which is also the global 
minimum. 

The state for 7^ is believed to be the quantum analog of the classical spiral phase. 
The fact that Neel order survives beyond the classically stable region is an example of the 
promotion of collinear order by quantum fluctuations, a phenomenon that has been observed 



in many other systems (see, e.g., Refs. 13lJ27l |). Thus, this collinear ordered state survives 



for the quantum case into a region where classically it is already unstable. Indeed, one can 
view this behavior more broadly as another example of the more general phenomenon of 
order by disorder— that we have briefly alluded to above, in which quantum fluctuations act 
to select and stabilize an appropriate type of order (that is typically collinear) in the face of 
classical degeneracy or near-degeneracy. 

It is also particularly interesting to note that the crossover from one minimum (0 = 0, 
Neel) solution to the other (0 7^ 0, spiral) appears (except for the LSUB2 case) to be 
quite abrupt (for all other LSUBn cases with even n > 2) at this point (and see Figs. 
|2] and [3]). Thus, for example, for the LSUB6 case the spiral pitch angle appears to 
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jump discontinuously from a zero value on the Neel side (k < 0.788) to a value of about 
O.IStt as the transition point into the spiral phase is crossed. This behavior is a clear 
first indication of a phase transition. Based on this evidence alone it would also appear 
that this transition is first-order, by contrast with the second-order nature of its classical 
counterpart. Such a situation where the quantum fluctuations change the nature of a phase 
transition qualitatively from a classical second-order type to a quantum first-order type 
has also been seen previously in the comparable spin-1/2 HAF model that interpolates 
continuously between the square and honeycomb lattices.— However, due to the extreme 
insensitivity of the energy to the pitch angle near the phase transition, as discussed above, 
we cannot rule out a continuous but very steep rise in pitch angle as the transition from the 
Neel phase into the spiral phase is transversed. All of the available evidence to date indicates 
that the transition at is subtle and may actually be second-order. Further evidence for 
the position and nature of the Neel-spiral quantum phase transition also comes from the 
behavior of the on-site magnetization that we discuss below. 

Before doing so, however, we wish to make some further observations on Figs. [2] and 
[3l We note first from Fig. [3] that in the case k = 1 (Ji = 1, = 1), corresponding to 
the spin-1/2 HAF on the triangular lattice, all of the CCM LSUBn approximations give 
precisely the classical value = f for the spiral angle, which corresponds to the correct 120° 
three-sublattice ordering. The fact that all LSUBn approximations give exactly this value 
is a consequence of us defining the LSUBn configurations on the triangular lattice (rather 
than the square lattice), and is a reflection of the exact triangular symmetry that is thereby 
preserved by our approximations. It is also interesting to note that for values of k > 1 the 
quantum spiral angle approaches the asymptotic (k oo) value of | much faster than 
does the classical angle. This is a first indication again, in this limit, of quantum fluctuations 
favoring collinear order (along the weakly coupled chains in this limit). 

We note from Fig. [2] that for certain values of J2 (or, equivalently, k) CCM solutions at 
a given LSUBn level of approximation (viz., LSUB6 in Fig. [2]) exist only for certain ranges 
of spiral angle 0. For example, for the pure square-lattice HAF (k = 0) the CCM LSUB6 
solution based on a spiral model state only exists for < < 0.147r. In this case, where 
the Neel solution is the stable ground state, if we attempt to move too far away from Neel 
coUinearity the CCM equations themselves become "unstable" and simply do not have a 
real solution. Similarly, we see from Fig. [2] that for k = 1.4 the CCM LSUB6 solution exists 
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FIG. 4: (Color online) Ground-state energy per spin versus J2 for the Neel and spiral phases of the 
spin-1/2 J1-J2 Hamiltonian of Eq. ([1]) with Ji = 1. The CCM results using the spiral model state 
are shown for various LSUBn approximations (n = {2, 4, 6, 8}) with the spiral angle <j) = (phSVBn 
that minimizes -ELSUBn (</>)• We also show the n — > 00 extrapolated result from using Eq. 

only for 0.247r < < O.Svr. In this case the stable ground state is a spiral phase, and now if 
we attempt to move too close to Neel collinearity the real solution terminates. 

Such terminations of CCM solutions are very common and are very well documented.— 
In all such cases a termination point always arises due to the solution of the CCM equations 
becoming complex at this point, beyond which there exist two branches of entirely unphysical 
complex conjugate solutions.— In the region where the solution reflecting the true physical 
solution is real there actually also exists another (unstable) real solution. However, only 
the (shown) upper branch of these two solutions reflects the true (stable) physicsl ground 
state, whereas the lower branch does not. The physical branch is usually easily identified 
in practice as the one which becomes exact in some known (e.g., perturbative) limit. This 
physical branch then meets the corresponding unphysical branch at some termination point 
(with infinite slope on Fig. [2]) beyond which no real solutions exist. The LSUBn termination 
points are themselves also reflections of the quantum phase transitions in the real system, 
and may be used to estimate the position of the phase boundary,™ although we do not do 
so for this first critical point since we have more accurate criteria discussed below. 

Thus, in Figs. H] and [5] we show the CCM results for the gs energy and gs on-site mag- 
netization, respectively, where the helical state has been used as the model state and the 
angle chosen as described above. For both quantities we show the raw LSUBn data for 



LSUB2 ■■ 
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n = {2, 4, 6, 8} and the extrapolated (LSUBcx)) results obtained from them by using Eqs. ([2]) 
and ([3]) respectively. Firstly, the gs energy (in Fig. H]) shows signs of a (weak) discontinuity 
in slope at the critical values discussed above. These values for themselves depend 
weakly on the approximation level. 

Secondly, the gs magnetic order parameter in Fig. O shows much stronger and much 
clearer evidence of a phase transition at the corresponding values previously observed 
in Fig. [3l The extrapolated value of M shows clearly its steep drop towards a value very 
close to zero at = 0.80 ± 0.01, which is hence our best estimate of the phase transition 
point. From the Neel side (k < the magnetization seems to approach continuously 
a value M = 0.025 ± 0.025, whereas from the spiral side {k > there appears to be 
a discontinuous jump in the magnetization The transition at k = k.-, thus 

appears to be (very) weakly first-order but we cannot exclude it being second-order, since 
we cannot rule out the possibility of a continuous but very steep drop to zero of the on-site 
magnetization from the spiral side of the transition, for the same reasons as 

enunciated above in connection with our discussion of Fig. [2] and [3l We find no evidence 
at all for any intermediate phase between the quasiclassical Neel and spiral phases. These 
results may be compared with those for the same model of Weihong et al.^ who used a 
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FIG. 5: (Color online) Ground-state magnetic order parameter (i.e., the on-site magnetization) 
versus J2 for the Neel and spiral phases of the spin-1/2 J1-J2 Hamiltonian of Eq. ([1]) with Ji = 1. 
The COM results using the spiral model state are shown for various LSUBn approximations (n = 
{2, 4, 6, 8}) with the spiral angle (j) = i;^LSUBn that minimizes £^LSUBn (</>) • We also show the n — > 00 
extrapolated result from using Eq. ([3]). 
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linked-cluster series expansion technique. They found that while a nonzero value of the Neel 
staggered magnetization exists for < k < 0.7, the region 0.7 < k < 0.9 has zero on-site 
magnetization, and for n > 0.9 they found evidence of spiral order. Nevertheless, their 
results came with relatively large errors, especially for the spiral phase, and we believe that 
our own results are probably intrinsically more accurate than theirs. 

As a further indication of the accuracy of our results we show in Table [IT] data for the 
two cases of the spin-1/2 HAF on the square lattice (/« = 0) and on the triangular lattice 
(k = 1). For both cases we present our CCM results in various LSUBn approximations 
(with 2 < n < 8) based on the triangular lattice geometry using the spiral model state, 
with = for the square lattice and = | for the triangular lattice. Results are given 
for the gs energy per spin E/N, and the magnetic order parameter M. We also display our 
extrapolated {n oo) results using the schemes of Eqs. ([2]) and ([3]) with the three data 
sets n = {2, 4, 6, 8}, n = {4, 6, 8} and n = {3, 5, 7}. The results are seen to be very robust 
and consistent. For comparison we also show the results obtained for the two lattices using 
quantum Monte Carlo (QMC) methods^*^ and linked-cluster series expansions.—"^ For the 
square lattice there is no dynamic frustration and the Marshall-Peierls sign rule^^ applies, so 
that the QMC "minus-sign problem" may be circumvented. In this case the QMC results^^ 
are extremely accurate, and indeed represent the best available for the spin-1/2 square-lattice 
HAF. Our own extrapolated results are in complete agreement with these QMC benchmark 
results, as found previously (see, e.g., Ref. j25|] and references cited therein), even though the 
LSUBn configurations are defined here on the triangular lattice geometry. Thus, we note that 
whereas the individual LSUBn results for the spin- 1 / 2 square-lattice HAF do not coincide 
with previous results for this model (see, e.g., Ref. [25|) because previous results have been 
based on defining the fundamental LSUBn configurations on a square-lattice geometry rather 
than on the triangular-lattice geometry used here, the corresponding LSUBoo extrapolations 
in the two geometries are in complete agreement with each other. 

By contrast, the nodal structure of the gs wave function is not exactly known for the 
spin-1/2 triangular-lattice HAF, and the QMC minus-sign problem cannot now be avoided 
for such frustrated spin systems. The QMC results showii^^ for the triangular lattice in 
Table [TTl were performed using a Green's function Monte Carlo method with a fixed-node 
approximation that was then relaxed in a controlled but approximate way using a stochastic 
reconfiguration technique. For the triangular lattice case we also show results in Table [III 
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TABLE II: Ground-state energy per spin and magnetic order parameter (i.e., the on-site magneti- 
zation) for the spin-1/2 HAF on the square and triangular lattices. We show COM results obtained 
for the Ji-J'2 model with Ji > 0, using the spiral model state in various LSUBn approximations 
defined on the triangular lattice geometry, for the two cases k = J2/J1 = (square lattice HAF, 
= 0) and K = 1 (triangular lattice HAF, (j) = We compare our extrapolated (n — > 00) results 
using Eqs. ^ and ([3]) and various sets of LSUBn data with other calculations. 



Method ■ 


E/N 


M 


E/N 


M 












square 


(k = 0) 


triangular (k = 1 ) 


LSUB2 


-0.64833 


0.4207 


-0.50290 


0.4289 


LSUB3 


-0.64931 


0.4182 


-0.51911 


0.4023 


LSUB4 


-0.66356 


0.3827 


-0.53427 


0.3637 


LSUB5 


-0.66345 


0.3827 


-0.53869 


0.3479 


LSUB6 


-0.66695 


0.3638 


-0.54290 


0.3280 


LSUB7 


-0.66696 


0.3635 


-0.54502 


0.3152 


LSUB8 


-0.66816 


0.3524 


-0.54679 


0.3018 


Extrapolations 


LSUBoo " 


-0.66978 


0.3148 


-0.55113 


0.2219 


LSUBoo ^ 


-0.66974 


0.3099 


-0.55244 


0.1893 


LSUBoo " 


-0.67045 


0.3048 


-0.55205 


0.2085 


QMC '^•'^ 


-0.669437(5) 0.3070(3) ■ 


-0.5458(1) 


0.205(10) 


SE 


-0.6693(1) 


0.307(1) ■ 


-0.5502(4) 


0.19(2) 



" Based on n = {2,4,6,8} 
* Based on n = {4, 6, 8} 
Based on n = {3, 5, 7} 

QMC (Quantum Monte Carlo) for square lattice^^ 
^ QMC for triangular lattice2S 
f SE (Series Expansion) for square lattice^ 
^ SE for triangular lattice^ 
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from a large-scale calculation using a linked-cluster series expansion.— For such frustrated 
systems this method, along with our CCM, is probably among the most accurate available. 
We see that in this case our results for the gs energy are in good agreement with the series 
expansion results, whereas the QMC estimate for the energy is almost certainly too high, 
and its quoted error hence erroneous. Our best estimate for the sublattice magnetization 
in this case, M = 0.20 ± 0.02, is in complete agreement with the best available by all other 
methods. 

The good agreement, both with respect to internal consistency checks using different 
extrapolations and with respect to other methods, for the gs properties of both the above 
models, gives us considerable confidence in our results for the spin- 1/2 Ji-Jg model for all 
values of K = Ji. 

We also comment on the k — »• oo (decoupled spin-1/2 ID HAF chains) limits of Figs. H] 
and [51 Firstly, Fig. H] shows that at large the extrapolated energy per spin approaches 
the value E/N = — 0.4431J2 which is the same as the exact result^ from the Bethe ansatz 
solution.— By contrast, the extrapolated magnetic order parameter at large seems to 
approach a constant value M ^ 0.10, by contrast with the exact value of zero^*^ in this 
limit. We note that the ID anisotropic XX Z chain with anisotropy parameter A has an 
essential singularity for M ^ at the isotropic point A — > and this is extremely difficult 
to mimic in any truncated numerical calculation. We note, however, that in the regime 
1 < K < 2 the order parameter M decreases almost linearly, and if this linear decrease were 
to be extended M would become zero at a value k, ~ 3.5. 

We turn finally to our CCM results based on the stripe state as CCM gs model state |$). 
The LSUBn configurations are again defined with respect to the triangular lattice geometry, 
exactly as before. Results for the gs energy and magnetic order parameter are shown in Figs. 
|6] and [7] respectively for the collinear stripe phase. They are the precise analogs of Figs. H] 
and [5] for the (Neel and) spiral phases. We see from Fig. [6] that some of the LSUBn solutions 
based on the stripe state show a clear termination point nt of the sort discussed previously, 
such that for k < no real solution for the stripe phase exists. In particular the LSUB6 
and LSUB8 solutions terminate at the values shown in Table IIIII As is often the case the 
LSUB2 solution does not terminate, while the LSUB4 solution shows a marked change in 
character around the value n ~ 0.880 that is not exactly a termination point (but, probably, 
rather reflects a crossing with another unphysical solution). In any event, the LSUB4 data 
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FIG. 6: (Color online) Ground-state energy per spin versus J2 for the stripe-ordered phase of the 
spin-1/2 J1-J2 Hamiltonian of Eq. ([T|) with Ji = 1. The CCM results using the stripe model 
state are shown for various LSUBn approximations (n = {2,4,6,8}). We also show the n — > 00 
extrapolated result from using Eq. 

are not shown below this value in Figs. [6] and [71 

The large k limit of the energy per spin results of Fig. [6] again agrees with the exact ID 
chain result of E/N = — 0.4431J2, just as in Fig. H] for the spiral phase. However, the most 
important observation is that for all LSUBn approximations with n > 2 the curves for the 
energy per spin of the stripe phase cross with the corresponding curves (i.e., for the same 

0.4 
0.35 
0.3 

0.25 

S 0.2 

0.15 
0.1 
0.05 


1 2 3 4 5 

FIG. 7: (Color online) Ground-state magnetic order parameter (i.e., the on-site magnetization) 
versus J2 for the stripe-ordered phase of the spin-1/2 J1-J2 Hamiltonian of Eq. ([T]) with Ji = 1. 
The CCM results using the stripe model state are shown for various LSUBn approximations (n = 
{2,4,6,8}). We also show the n — > 00 extrapolated result from using Eq. 
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TABLE III: The parameters Kg (the crossing point of the energy curves for the stripe and spiral 
phases) and Kt (the termination point of the stripe state sohition) in various LSUBn approximations 
defined on the triangular lattice geometry, for the spin-1/2 J1-J2 model, with n = J2/J1, Ji > 0. 
The "LSUBoo" extrapolations are explained in the text. 

LSUBn 





Kg 




LSUB2 


00 




LSUB4 


4.555 


(0.880) 


LSUB6 


3.593 


0.970 


LSUB8 


3.125 


1.150 


TSUBoo" 


L69 ± 0.03 


1.69 



value of n) for the energy per spin of the spiral phase at a value that we denote as He- Thus, 
for K, < He the spiral phase is predicted to be the stable phase (i.e., lies lowest in energy), 
whereas for k > Kg the stripe phase is predicted to be the stable ground state. We thus 
have a clear first indication of another (first-order) quantum phase transition in the spin-1/2 
J1-J2 model at a value n = k^^- Figure [8] shows the energy difference between the stripe and 
spiral states for various LSUBn calculations, and some indicative values near the crossing 
values also shown in Table UTIl 

As we remarked in Sec. [TTl the stripe phase is never the stable classical ground state 
since it always lies higher in energy than the spiral phase. Indeed, it is easy to show that 
the classical gs energy per spin, written as e = E{Ns^J'2)~^-, has the value e^f''"' = — 1 
(independent of Ji) for the stripe phase, and e^f^'^'^' = — 1 — for the spiral phase in 
the regime > | where the spiral phase exists classically. Thus, the difference in energy 
between the two phases classically is e^^^"^^^ — g^t"P<= = /^g^^ = < 0- By contrast, what 

we have shown at the quantum level is that quantum fluctuations can change the sign of Ae 
at some critical value k = such that the collinear stripe phase becomes stabilized for all 
K > Kc2- However, the energy differences Ae in this regime are found to be extremely small, 
as may be seen from Fig. [HI and hence the stripe phase is predicted to be very fragile against 
small perturbations or thermal fluctuations, for example. Nevertheless, we should stress 
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that the LSUBn energy differences displayed in Fig. [8] are well within (by several orders of 
magnitude) the margins of error in our individual calculations. 

Clearly, the LSUBn energy crossing points Ke(n) at which AcLsuBn = provide a measure 
of Kcj. In the past we have found that a simple linear extrapolation, k = cq + Cin~^, yields 
a good fit to such critical points, and this seems to be the case here too. The corresponding 
"LSUBcxd" estimate from the LSUBn data of Table UTTl with n = {4, 6, 8} gives an estimate 
Kc'z ~ 1-69 ± 0.03, where the error is the standard deviation in the fit. A similar linear 
extrapolation on the LSUBn stripe-phase termination points with n = {6, 8} gives a 
second estimate k,c2 ~ 1-69, in remarkably good agreement with the first estimate. We note 
too that we can, of course, also obtain another estimate of Kg from the crossing point of 
the two extrapolated (LSUBcxo) gs energy curves for the spiral and stripe-ordered phases 
(i.e., using Eq. ([2]) in each case), shown in Figs. [Hand [6] respectively, and it is this LSUBoo 
result that is displayed in Fig. [HI Since the two curves cross at such a very shallow angle, 
and since both curves are extrapolated completely independently of one another, we expect 
that this estimate for Ke is perhaps intrinsically less accurate than the one discussed above. 
Nevertheless, it is very gratifying that the value so obtained for the crossing point of the two 
extrapolated (LSUBoo) gs energy curves, namely ^ 2.19, is rather close to the previous 

0.0002 
0.0001 


< -0.0001 
-0.0002 
-0.0003 
-0.0004 

1 1.5 2 2.5 3 3.5 4 4.5 5 

h' 

FIG. 8: (Color online) Difference between the ground-state energies per spin (e = E/N) of the 
spiral and stripe phases (Ae = e^^"^^ — e^*'''P*') versus J2 for the spin-1/2 J1-J2 Hamiltonian of 
Eq. ([1]) with Ji = 1. The CCM results for the energy difference using the stripe and spiral model 
states for various LSUBn approximations (n= {4, 6, 8}) are shown. We also show the n ^ 00 
extrapolated result from using Eq. ([2]) for the two phases separately. 
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TABLE IV: Comparision of the gs energy per spin of the stripe and spiral phases of the spin-1/2 
Ji-J'2 model in various LSUBn approximations (n = {2, 4, 6, 8}) at some illustrative values of J2 
(with Ji = 1). 

E/N 

LSUBn JL 







Li IJJC 


QTMTa 1 
0^11 Chi. 


9 


1 n 




-4:. LXjOOLO 


4 


4.4 


-1.927326 


-1.927337 


4 


4.6 


-2.014224 


-2.014219 


4 


4.8 


-2.101154 


-2.101136 


6 


3.4 


-1.508386 


-1.508406 


6 


3.6 


-1.595632 


-1.595632 


6 


3.8 


-1.682968 


-1.682952 


8 


2.9 


-1.295557 


-1.295584 


8 


3.1 


-1.382665 


-1.382667 


8 


3.3 


-1.469937 


-1.469924 



value, considering that the two curves are almost parallel to one another in the crossing 
region. The difference in the two estimates for Kg of about 1.7 and 2.2 is itself an indication 
of the error in these estimates in this incredibly difficult regime where the two phases lie so 
close in energy to one another. Nevertheless, we reiterate that the CCM retains sufficient 
accuracy for both the spiral and stripe phases individually (both in the raw LSUBri data 
and in the extrapolated LSUBcxd results for each phase) to ensure that, even though the 
corresponding estimates for the energy difference Ae are small in the crossing regime and 
for values of k > Kc, they are still sufficiently large to be well within our limits of accuracy. 

Finally, from Fig. [7] we see the corresponding results for the magnetic order parameter of 
the stripe phase. The large-K limit of uncoupled ID spin-1/2 chains is identical to that for 
the spiral phase, and hence suffers from the same (known) problem of not giving the correct 
M = result in this limit. However at smaller values of k the order parameter decreases 
and becomes zero at some critical value Km{n) that depends on the LSUBn approximation. 
The extrapolated [n —>■ 00) result obtained in the usual way from Eq. ([3]) is also shown in 
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Fig. [7] and is seen to become zero at a value Km ~ 1.43. Since the phase transition at k = k,c2 
is clearly of first-order type from the energy data, the magnetization data provide us only 
with an inequality k,c2 ^ '^m- 

In summary, although it is difficult to put firm error bars on our results for our predicted 
second critical point, our best current estimate, based on all the above results, is = 
1.8 ±0.4. 

V. DISCUSSION AND CONCLUSIONS 

In this paper we have used the CCM to study the infiuence of quantum fluctuations on 
the zero-temperature gs phase diagram of a spin-half Heisenberg antiferromagnet, the J1-J2 
model, defined on an anisotropic 2D lattice. We have studied the case where the NN Ji 
bonds are antiferromagnetic (Ji > 0) and the competing = nJi bonds have a strength k 
that varies from k = (corresponding to the spin-half HAF on the square lattice) to k — ^ cx3 
(corresponding to a set of decoupled spin-half ID HAF chains), with the spin-half HAF on 
the triangular lattice as the special case k = 1 in between the two extremes. 

Whereas at the classical level the model has only two stable gs phases, one with Neel order 
for K < = 0.5 and another with spiral order for k, > k^i, the quantum fiuctuations for 
the spin-half model can stabilize a third non-classical phase with coUinear stripe ordering at 
sufficiently high values of k. Thus, we find two quantum phase transitions, both seemingly 
first-order. The first at k,c\ = 0.80 ± 0.01 separates a phase with classical Neel ordering 
for K < from a phase with helical ordering for n > k.^- This latter phase includes the 
case K = 1 of the spin-half triangular-lattice HAF with the standard 120° three-sublattice 
quasiclassical ordering. By contrast with the classical second-order transition at = 0.5, 
the quantum phase transition at ~ 0.80 appears to be weakly first-order in nature, 
although we cannot exclude it from being second-order. The Neel order thus survives into 
a region 0.5 < k < where it is classically unstable. This is an example, among many 
others, of the widely observed phenomenon that quantum fiuctuations tend to favor collinear 
ordering. 

A second quantum phase transition is predicted by our calculations to occur at k,c2 = 
1.8 ± 0.4. It separates the helical phase for k, < k,c2 from a phase with collinear stripe 
ordering at k > k,c2- However in this latter region, k > Kcj, the stripe and spiral phases are 
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extremely close in energy, and hence the stripe phase may be very fragile against external 
perturbations or thermal fluctuations, for example. The existence of the coUinear stripe 
phase seems to rely again on the fact that quantum fluctuations favor collinear ordering. 
An alternative, but essentially equivalent, explanation starts by looking at the large-K limit 
of uncoupled ID HAF chains, for which all states with spins on different chains orienting 
randomly with respect to each other are degenerate in energy in the k — oo limit. Then, 
as the interchain coupling Ji = J^k is slowly increased from zero, what seems to occur 
is that the relative orientations are locked into collinearity by the familiar phenomenon of 
order by disorder.— As a somewhat technical aside at this point we note, however, that while 
the present J1-J2 niodel attains (collinear) order from an otherwise disordered set of spin 
chains, the order by disorder phenomenon seemingly responsible here differs in an important 
respect from its archetypal realization in the J1-J2 model. In this latter case each of the 
inter-penetrating sublattices is characterized by a classical magnetization vector, and the 
classical disorder emanates from the independence of the energy on the mutual orientation 
of the magnetization vectors for the two sublattices. By contrast, for the present J1-J2 
model, the corresponding phenomenon is more subtle since the spin chains that now order 
(collinearly) are themselves quantum-critical. Hence, what is established now is not only 
that the relative orientations of the magnetizations on different chains become locked into 
collinear order, but also the very fact that such a classical notion is actually appropriate. 

We note that our own calculations are weakest precisely in the k — 00 limit where, 
although our result for the gs energy of the chains is excellent, the subtle quantum-critical 
ordering of the Bethe ansatz solution that results in the value M = for the staggered 
magnetization,^ is not exactly reproduced. Nevertheless, with the single exception of our 
inability to reproduce exactly this very singular and non-analytic result for the on-site 
magnetization,— the CCM calculations are remarkably robust and accurate over the rest 
of the parameter space. Thus, our results for the gs energy and on-site magnetization have 
provided a set of independent checks that lead us to believe that we have a self-consistent 
and coherent description of this interesting model. 

Before concluding we compare our results with those from previous calculations on the 
same model. The simplest such studies have utilized lowest-order (or linear) spin-wave 
theory (LSWT)^'^ and Schwinger boson mean-field theory (SBMFT).^ The independent 
LSWT calculations of Trumper- and Merino et ai— both found a continuous second-order 
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phase transition from a Neel-ordered to a spiral-ordered phase at precisely the classical 
value = 0.5, at which the magnetization order parameter approaches zero continuously 
from both sides. Although LSWT is known to give a reasonable description of the spin-1/2 
Heisenberg antiferromagnet on both the square lattice (k = 0) and the triangular lattice 
(k = 1), it is clearly unable to model the intermediate regime accurately. This is particularly 
true around the point of maximum classical frustration at k, = 1/2, for which one fully 
expects, as mentioned previously, that Neel order is preserved to higher values of k than 
pertain classically, as found both by earlier more accurate studies, such as those using series 
expansion (SE) techniques,- and by us in the present paper. 

Similar shortcomings of spin-wave theory (SWT) have been noted by Igarashi^ in the 
context of the related spin-1/2 J1-J2 model on the square lattice, discussed briefly in Sec. 
[H He showed that whereas its lowest-order version (LSWT) works well when J2 = 0, it 
consistently overestimates the quantum fluctuations as the frustration J2/J1 increases. In 
particular he showed by going to higher orders in SWT in powers of 1/s where s is the spin 
quantum number and LSWT is the leading order, that the expansion converges reasonably 
well for J2/J1 ^ 0.35, but for larger values of J2/J1, including the point J2/J1 = 0.5 of 
maximum classical frustration, the series loses stability. He also showed that the higher- 
order corrections to LSWT for J2/J1 ^ 0.4 make the Neel-ordered phase more stable than 
predicted by LSWT. He concluded that any predictions from SWT for the spin-1/2 J1-J2 
model on the square lattice are likely to be unreliable for values J2/J1 ^ 0.4. It is likely that 
a similar analysis of the SWT results for the present spin-1/2 J1-J2 niodel on the square 
lattice would reveal similar shortcomings of LSWT as the frustration parameter k, = J'2/ Ji 
is increased. 

By contrast with the above LSWT results for the spin-1/2 Ji-J'2 model on the square 
lattice, a SBMFT analysis^ shows a continuous transition from a collinear Neel phase to a 
spiral phase at a value ^ 0.6, but with a non- varnishing magnetization, M ^ 0.175, 
at the critical point. The discrepancy with our own results (viz., ^ 0.80 with either a 
vanishing or very small magnetization, M ^ 0.025 ± 0.025, at the critical point) is almost 
certainly again due to the lowest-order nature of the mean-fleld approach, and particularly 
the complete neglect at the SBMFT level of Gaussian fluctuations. 

Turning to the second phase transition we note that LSWT predicts that the magneti- 
zation in the spiral phase of the present spin-1/2 Ji-J'2 model vanishes at a value k ~ 3.70. 
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The authors of these results^*^ took that to indicate the possible existence of a disordered 
phase for k > 3.70. Even ignoring the probably unreliable nature of LSWT results for such 
high values of k, as noted above, the safer conclusion based on them is that the spiral phase 
simply becomes unstable for k > 3.70. Indeed, SWT results are always based on a particu- 
lar choice of phase, usually based on a classical (s — oo) analog, just as are our own CCM 
calculations. In both cases the vanishing of an order parameter only signals a phase transi- 
tion to another state, but in the absence of another calculation based on another (classical) 
hypothesized phase, no conclusion about the adjacent phase can be drawn. It is precisely 
such an additional calculation based on a stripe-ordered phase that has been done in the 
present work, which has shown the onset of this phase (at a lower-energy than the spiral 
phase) for k > /tcj ~ 1.8 ±0.4. No such calcultions were attempted in either the LSWT-i^ or 
the SBMFT- cases, and hence no direct comparison can be drawn with own results for this 
second phase transition, except to say that the LSWT result with its expected uncertainties 
discussed above is not incompatible with our own conclusions. 

As has been noted elsewhere,- high-order CCM results of the sort presented here are 
believed to be among the best available for such highly frustrated spin-lattice models. Many 
previous applications of the CCM to unfrustrated spin models have given excellent quan- 
titative agreement with other numerical methods (including exact diagonalization (ED) of 
small lattices, quantum Monte Carlo (QMC), and series expansion techniques). A typical 
example is the spin-half HAF on the square lattice, which is the k = limit of the present 
model (and see Table HTl) . It is interesting to compare for this k = case, where comparison 
can be made with QMC results, the present CCM extrapolations of the LSUBn data for the 
infinite lattice to the n — oo limit and the corresponding QMC or ED extrapolations for the 
results obtained for finite lattices containing spins that have to be carried out to give the 
N oo limit. Thus, for the spin-1/2 HAF on the square lattice the "distance" between the 
CCM results for the ground-state energy per spin^^ at the LSUB8 (LSUBIO) level and the 
extrapolated LSUBcxd value is approximately the same as the distance of the corresponding 
QMC result^ for a lattice of size A^ = 12 x 12 (A^ = 16 x 16) from its A^ cx) limit. The 
corresponding comparison for the magnetic order parameter M is even more striking. Thus 
even the CCM LSUB6 result for M is closer to the LSUBoo limit than any of the QMC 
results for M for lattices of A^ spins are to their N (yo limit for all lattices up to size 
A^ = 16 X 16, the largest for which calculations were undertaken.— Such comparisons show. 
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for example, that even though the "distance" between our LSUBn data points for M and 
the extrapolated {n oo) LSUBoo result shown in Fig. [5] may, at first sight, appear to 
be large, they are completely comparable to or smaller than those in alternative methods 
(where they can be applied). Furthermore, where such alternative methods can be applied, 
as for the spin-1/2 HAF on the square lattice, the CCM results are in complete agreement 
with them. 

By contrast, for frustrated spin- lattice models in two dimensions both the QMC and ED 
techniques face formidable difficulties. These arise in the former case due to the "minus-sign 
problem" present for frustrated systems when the nodal structure of the gs wave function is 
unknown, and in the latter case due to the practical restriction to relatively small lattices 
imposed by computational limits. The latter problem is exacerbated for incommensurate 
phases, and is compounded due to the large (and essentially uncontrolled) variation of the 
results with respect to the different possible shapes of clusters of a given size. 

Thus, for highly frustrated spin-lattice models like the present J1-J2 model, the best 
alternative numerical method to the CCM is the linked- cluster series expansion (SE) 
technique .-'^'^'^'^ The SE technique has also been applied to the present model.-'^ The 
earlier study^ mainly dealt with the Neel and spiral phases. Unlike in that work we find 
no evidence at all for an intermediate (dimerized) phase between the Neel and spiral phases 
in the parameter regime 0.7 ^ k ^ 0.9. The very recent SE study^^ was motivated by the 
prediction of Starykh and Balents^ for the existence of a stable coUinear stripe-ordered gs 
phase for values of n above some critical value that they did not calculate. The SE study 
showed that although the collinear stripe phase was stabilized for large values of k relative to 
the classical result, nevertheless in their calculations the non-collinear helical phase was still 
always lower in energy. Hence, they could not confirm the existence of the stripe-ordered 
phase. They concluded by suggesting that further unbiased ways of studying the competi- 
tion between the spiral and stripe phases would be useful. We believe that the present CCM 
calculations provide exactly such unbiased results, which now do indeed appear to confirm 
the prediction of Starykh and Balents. 

We end by remarking that it would also be of interest to repeat the present study for 
the case of the spin-one J1-J2 model. The calculations for this case are more demanding 
due to an increase at a given LSUBn level of approximation in the number of fundamental 
configurations retained in the CCM correlation operators. Nevertheless, we hope to be able 
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to report results for this system in the future. 
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